the velocities still depend only on the index k of the hit boundary simplex, but where 
they are not restricted by this constraint. By the Perron-Frobenius theorem there also 
exists an unique probability measure but its algebraic expression is not as simple as in 
the case above, especially for larger values of n. 

A further class of dynamical systems is obtained by choosing the velocities not only 
dependent on the faces at the actual "reflection" as in the present work, but also 
depending on the particular position where the trajectory hits the boundary [8]. 
In such systems one also observes simplex-like invariant sets with singularities but 
their dynamical properties are quite different. In [8] we will introduce briefly the 
generalized dynamical systems and will also discuss possible applications to chemical, 
bio-technological manufacturing and socio-economical systems. 
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Theorem 3.3 Let G : dS n — > dS n be defined by equation (4). Let also the dynamics 
be determined by the point p. Then the entropy h^ corresponding to the map G is 

1 n 1 — 
hfj, = - pi(l-pi) log with (24) 

n 

d = pi t 1 ~ ^')- 

8 = 1 

Proof: Let us consider the finite partition V = {Si jfl -i \ i = 1, . . . , n} of dS n . Since the 
boundaries dSi >n -i of its elements are invariant w.r.t. the map G, i.e. the boundaries 
of the subsets Si >n -i can not be mapped in interiors, the partition V is a Markov 
partition[5]. Markov partitions for which all intersections Si >n -i r\G(Sj >n -i) and Si >n -i D 
G~ x (S hn -\) are connected are generating partitions [6]. This is indeed the case for the 
above system. Also our partition has the generating property. The symbolic dynamics 
a (c.f. sec. 2) is thus equivalent to a discrete time Markov process, i.e. the probability 
of hitting face S{ ;a -\ at time t m is only dependent on the preceding face hit at t m _i. 
For the complete statistical characterization of the system, let pi = pi (1 — pi)/d, i = 
1, . . . , n, be the weight of the z-th subsimplex. The (conditional) transition probabilities 
for transitions j ' — > i are given by the matrix pij = (dj)i. The shift a on the space of 
sequences is a measure preserving transformation for the Markov measure defined by 
Pij and initial vector pj, satisfying pi = ^2jPij Pj- Also the transformation a is ergodic 
because the chain is irreducible, i.e. for all pairs of states there exists a m > 
with p\f > (e.g. m = 2). Since the entropy of an ergodic Markov shift is [6] 

hf, = ~ Y V] Pij log Pij, (25) 

i,j 

by inserting the above measures and doing some algebraic simplifications we get equa- 
tion (24). 

Hence for n > 2 we get positive entropy. Since all transitions i — > j with i ^ j are 
allowed, while z — > z is forbidden, the topological entropy is h top = log(n — 1). 

4 Summary and outlook 

Inspired by hybrid dynamical systems recently investigated in modelling chemical man- 
ufacturing, we introduced dynamical systems on n-simplex geometries, generalizing 
thereby the results obtained in [3] to arbitrary n. 

By use of generating partitions it was possible to define symbolic Markov dynamics. 
Due to the ergodicity of the corresponding Markov shift a it was possible to compute 
the metric entropy with help of the underlying SRB measure. 

In the considered system the velocities = — p of the particle are not independent, 
but Vk — vi = ek — e\. One might ask what happens in the more general case where 
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Proof: The simplexes Sk, n -i an d S{ n _ 1 have in common an entire (n — 3)-dimensional 
subsimplex, namely the face opposing the corner dk of S\ n _ 1 . We call this face Sk-i, n -2- 
Therefore we have 



HS 1 ^) dist(4,^,.- 2/ 



MS k ,n-i) dist(e,,5' fc; ,,„_ 2 ) 



where 



dist(x, A) = min { \x — y | : y £ A}. (20) 



The vectors dk — y\ and e; — y 2 (where y\ and y 2 minimize the denumerator and 
denominator in eq.(19), respectively) are both in the linear space spanned by Sk, n -i 
and perpendicular to Sk-i, n -2- They are thus parallel vectors, and the ratio of their 
lengths is equal to the ratio of any one of their components. Taking in particular the 
/-th component (for which (y\)i = (1)2)1 = and (e;); = 1) we find eq.(18). 

With lemma 3.1 and eq.(17) we get immediately the following theorem: 

Theorem 3.1 Any measure A which has constant density fik,n-i on each subsimplex 
Sk, n -i evolves under the application of the map G according to 



^ [dk)l (ik,n-l- (21) 



k=l 



Note that by using ^™ =1 Pi = 1 and the fact that p is in the interior of S n , one easily 
verifies the identity ^™ =1 (^fc)/ = 1 for / = 1, . . . , n, proving that the matrix pik = (dk)i 
is indeed a stochastic matrix. 

The following theorem yields a fixed point of eq.(21) 
Theorem 3.2 The piece-wise constant measure 

X*(Si tn - 1 ) = -pi(l-pi), i = l,...,n, (22) 

a 

n 

d = E P^ l -Pi), (23) 

8 = 1 

is a properly normalized probability measure and invariant under G. 
Proof: By straightforward computation. 

For n = 3, this has been derived already in [3]. Indeed the above measure is absolutely 
continuous and also unique due to the Perron-Frobenius theorem [7]. Such measures 
are called SRB measures. Its corresponding density is given by the relations (15,16). 
Finally, we state our main result: 



6 



subsets Si >n -i. They are disjoint, apart from their edges which have measure zero. The 
measure of subsimplex i can be expressed as 

KSi,n-i) = I //(*) dA{x) (11) 

The probability density function ji(x) is called to be invariant w.r.t. the map G, if 

\(G-\B)) = \(B) (12) 

for all B C dS n . Suppose that the invariant probability density function // cannot be 
written as the superposition // = + (1 — a)fi 2 with < a < 1, /Ji / (J^- Then // is 
called indecomposable or ergodic. Regardless of the initial probability density function, 
provided it was absolute continuous w.r.t. Lebesgue, the system will asymptotically 
be described by the invariant probability density function fj* . 

Now let us construct the Frobenius-Perron operator and determine the invariant prob- 
ability density function fj* on dS n . By 1a(%) we denote the characteristic (or indicator) 
function for a set A simply defined by 

Then we consider the ansatz of a piece-wise constant density 

n 

M^) = ^2^,n-l ISi^ix), //,>-! G R + . (14) 

8 = 1 

For the particular sets Sk, n -i, k = 1, . . . , n, we find by the above definition 

MSk,n-l) = Hk,n-1 MSk,n-l) (15) 

and 

KSi, n -l) = ^lA(t-l)- ( 16 ) 

Here, A(A) denotes the Lebesgue measure of set A in dimension (n — 2) and S l kn _ 1 
is the subset of Sk, n -i which is mapped under G to Si >n -i, see fig.(l). Note that G~ x 
is not a single valued map. Invariance of the measure A requires then a fixed point 
equation for the densities: 

/V-i = s> T7^ \ / = l,...,n. (17) 

Lemma 3.1 In the above constructed state space and under the mentioned assumptions 
it is 

A 1 < ~'h,n-\, 



^^r = (d k ,e,). (18) 
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This gives (p^ — 1) a + 1 = 0, and thus 

1 



a = - , i.e. (7) 

1 - Pk 

d k = -. (p 1 ,p 2 ,...,0 k ,...,p n ). (8) 

1 - Pk 

We have thus proven the following lemma: 

Lemma 2.1 Let dk £ Sk, n -i be defined by eq.(8). It follows that 

dk 1 — > 4, (9) 

and the time step associated with this map is At = a, where a is given by eq.(7). 

We note that the trajectories starting with a point x £ Sk, n -i, see ec [-(4) are always 
parallel to the line between dk and through the point p } because of the same direc- 
tional vector of the trajectory. 

Finally let us define a symbolic dynamics. Each time t m when a face S Wrn>n -i is hit by 
the particle we write down the index u m of the actual face. Then any point x which 
gives rise to a trajectory with x(t ) = x and x(t m ) £ S Wrn>n -i is coded by an infinite 
sequence 

u = {u ,ui, . . .}, uj k £ {1, . . . ,n}. (10) 

The dynamics on the space of sequences is indeed defined by the the ordinary shift 
transformation u' = cr(u;), where u' k = LOk+i- For the computation of the entropy of G 
one can equivalently determine the entropy of the shift a if the partition is generating [6] . 

In the following section we first determine an invariant ergodic measure of the map G. 
Further we show that the partition beyond the transformation a is generating and that 
a represents an ergodic Markov shift. The entropy of the map G is then straightforward 
determined by computing the entropy corresponding to the Markov measure. 



3 Determining invariant measures and entropies 

For what follows we give here a short repetition of the theoretical framework about 
measures, probabilities and ergodicity. We consider the invariant set dS n . For each 
measurable subset B of dS n , the symbol A(B) = J B dA(x) denotes the Lebesgue mea- 
sure of B. We look for an invariant probability measure which is absolutely continuous 
with respect to A, and for its density p. In nonlinear dynamics a A-density p is called 
an invariant probability density function w.r.t. the dynamics of the system. In view 
of the natural given partition V = {Si jfl -i \ i = 1, . . . , n} of dS n we consider the special 
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and x(t m+ i) is determined unique from x(t m ). Precisely, suppose that the trajectory 
hits at time t m the fe-th face i.e. x(t m ) G Sk, n -ii then At m is determined by simple 
geometrical arguments and is 



At r 



mm 




(3) 



#k \ pi 

Then we formally write the Poincare map by 

x(t m+1 ) = G(x(t m )) = x(t m ) + (4 - p) At m , (4) 
where is corresponding to the face Sk, n -i stroken at time t m . 



X ▼ 



Figure 1: Illustration of the state space of a "strange" billiard inside a 3-simplex. The spreading of 
nearby trajectories is presented by the grey area. The transition point p parameterizes the dynamics 
G. 

Of special interest are the points (4 G Sk, n -i from which the trajectory hits the opposite 
corner. It is defined by 1 



4 = d k + (4 - p)a. 

(4 , 4) = o. 



a G R, and 



(5) 
(6) 



1 ( , } denotes the canonical scalar product in R n 
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are several relations to hybrid dynamical systems recently investigated in modelling 
chemical manufacturing [3] [4]. 

In section 3 we determine the invariant ergodic measure (also called SRB measure) 
of the dynamics for any dimension. A generating partition is considered and we find 
that its corresponding symbolic dynamics is ergodic. Further we determine the (KS or 
"metric") entropy and find that it is positive for all cases n > 2. Finally the topological 
entropy is derived. 

2 Dynamics of "strange" billiards 

As mentioned in the introduction, "strange" billiards are defined as special standard 
billiards where the reflection rule is substituted by a "strange" rule. For the definition, 
let S n be the standard n-simplex embedded in R n 

^Xi = l, x 3 >0, for j = l,...,)i . (1) 

i = l ' 

Further we denote by Sk, n -i the k-th face of the boundary of S n . They are (n — 1)- 
simplexes 

S k , n -i = {x G R n \ x k = 0, x G S n } , for k = l,...,n. (2) 

n 

Note that the boundary of S n is just the union dS n = [J Si >n -i of the faces. 

8 = 1 

As long as x is not on the boundary of S n , the dynamics is as in the ordinary billiard: 
(x,v) — (x + vt,v). But in contrast the reflection at the boundary is not specular. 
Instead, for each face Sk, n -i there is a fixed velocity v k under which the trajectory 
leaves the face, independently of the direction of the incident velocity. Furthermore, 
we assume that these velocities are not independent, in particular v k = e k — p where 
ek is the k-th canonical unit vector in R n and p an element of the interior of S n . Note 
that the dynamics is completely defined by the choice of p. 

A realization of this model is a set of n containers served by a single server. The 
content X{ in the z-th container decreases by a rate pi } while the server can refill with 

n 

unit rate. To start with, X{ = 1 and the server is in position i . It remains there 

8 = 1 

filling container i until another container becomes empty. At this time the server 
switches instantaneously and remains in the new position until the next container is 
emptied, and so on[3]. 

For the definition of a symbolic dynamics we first consider "Poincare" sections by 
restriction of the trajectory to the boundary of S n . Therefore let us denote by t m 
the time where the trajectory xit) strikes the boundary dS n at the ra-th times. The 
time interval At m = t m+ \ — t m between two successive hits of the boundary at x(t m ) 
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Abstract 

In the present work we investigate a new type of billiards defined inside of 
n-simplex regions. We determine an invariant ergodic (SRB) measure of the 
dynamics for any dimension. In using symbolic dynamics, the (KS or metric) 
entropy is computed and we find that the system is chaotic for all cases n > 2. 

1 Introduction 

Let Q be a bounded region in R n with piece-wise smooth boundary (n > 2). A billiard 
in Q is a dynamical system generated by the uniform linear motion of a material 
point inside Q with a constant velocity, and with reflection at the boundary such that 
the tangential component of the velocity remains constant and the normal component 
changes sign[l]. The phase space of a billiard consists of all possible pairs (q,v) where 
q G Q and v is the velocity vector, i.e. an element of the unit sphere in n-dimensional 

space. Let M be the set of the points x = (q,v), and F t : M > M the flow on M. 

This dynamics is most conveniently replaced by a map G t : dM — > dM which maps 
(q n G dQ } v n ) onto (q n +i G dQ } v n+ i). Here, q n is on the boundary of Q and v n is the 
velocity before the point hits the boundary, while v n+ i is the velocity after reflection 
and q n +i the next point where the boundary is hit. 

Several statements are proved about billiards moving in various geometries [1] [2]. For 
example, the entropy of billiards inside n-dimensional polyhedrons is zero, while some 
other billiards are among the simplest hamiltonian systems proven to be ergodic and 
mixing. Thus the investigations of billiards are interesting from the dynamical system 
point of view. Also, a number of other problems of physics can be reduced to the study 
of billiards. 

In the present work we investigate a new type of dynamical systems, so-called "strange" 
billiards, moving in n-simplex regions. They are similar to the billiards discussed above, 
but the reflection rule at the boundary is substituted by a "strange" rule defined in 
the following section. Although the dynamics will seem somewhat artificial, there 



